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ABSTRACT 

Herschel-PACS measurements of the rotational R(0) and R(l) HD lines in the atmospheres of Uranus and Neptune are analyzed 
in order to derive a D/H ratio with improved precision for both planets. The derivation of the D/H ratio includes also previous 
measurements of the R(2) line by the Short Wavelength Spectrometer on board the Infrared Space Observatory (ISO). The available 
spectroscopic line information of the three rotational transitions is discussed and applied in the radiative transfer calculations. The best 
simultaneous fit of all three lines requires only a minor departure from the Spitzer temperature profile of Uranus and a departure limited 
to 2K from the Voyager temperature profile of Neptune (both around the tropopause). The resulting and remarkably similar D/H ratios 
for Uranus and Neptune are found to be (4.4+0.4)xl0~^ and (4.1+0.4)xl0~^ respectively. Although the deuterium enrichment in 
both atmospheres compared to the protosolar value is confirmed, it is found to be lower compared to previous analysis. Using the 
interior models of Podolak et al. (1995), Helled et al. (2011) and Nettelmann et al. (2013), and assuming that complete mixing of 
the atmosphere and interior occured during the planets history, we derive a D/H in protoplanetary ices between (5.75-7.0)xlO~^ for 
Uranus and between (5.1-7.7)xl0~^ for Neptune. Conversely, adopting a cometary D/H for the protoplanetary ices between (15- 
30)xlO~^, we constrain the interior models of both planets to have an ice mass fraction of 14-32%, i.e. that the two planets are 
rock-dominated. 

Key words. Planets and satellites: individual: Uranus - Planets and satellites: individual: Neptune - Planets and satellites: atmospheres 
- Planets and satellites: interiors 



1. Introduction 

Among the light nuclides synthesized during the early evolution 
of the universe, deuterium is unique in its sensitivity to deter- 
mine the cosmological density of baryons. As early as in the 
pre-main sequence of stars, deuterium was burned up to ^He. 
The gas, which returns into the interstellar matter by stellar out- 
flows and supernova explosions, is then free of deuterium. On 
the other hand, no process is known to produce deuterium, there- 
fore its abundance is decreasing with time, particularly with pro- 
gressing star formation. The deuterium abundance as measured 
today thus provides a lower limit for its corresponding primor- 
dial value. The determination of accurate D/H ratios in the at- 
mospheres of the Giant Planets has therefore been a longstand- 
ing target of research, because their values enable constraining 
the D/H ratio in the part of our Galaxy where the Solar sys- 
tem formed. Moreover, the D/H ratio is known to increase in 



* Herschei is an ESA space observatory with science instruments pro- 
vided by European-led Principal Investigator consortia and with impor- 
tant participation from NASA. 



icy grains with decreasing temper ature d ue to ion-molec ule and 
grain-surface interactions (Watson [l974[ Brown & Millar IT9891) . 
Measuring the D/H ratio as a function of heliocentric distance in 
the Solar system therefore enables probing the temperature of 
formation of icy grains in the protoplanetary disk (Owen et al. 
1999; Hersant et al. 2001; Gautier & Hersant 2005 ), since the 
D/H value measured in the atmospheric gas can be linked to the 
D/H value in the protoplanetary ices as will be shown in sec- 
tion |5] 

The D/H ratio in the hydrogen of the atmospheres of Jupiter 
and Saturn is believed to be very close to the protosolar value, 
because the mass of their cores is negligible with respect to their 
total mass and because H2 is by far the main constituent of their 
atmospheres. Consequently, deuterium enrichment of the hydro- 
gen reservoir through ices played only a minor role during their 
formation. However, the atmospheres of Uranus and Neptune are 
expected to have atmospheres enriched in deuterium. According 
to Guillot ( 1999 ), their cores (which in the models are composed 
of 25% rock and 60-70% of ice) make up for more than half of 
the total mass and mixing of deuterium-enriched icy grains and 
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planetesimals with the hydrogen envelope during their forma- 
tion (Hubbard & McFarlane 11980b must have led to a substan- 
tially larger D/H ratio in their atmospheres as compared to the 
protosolar value. 

Molecular hydrogen represents the major fraction of the at- 
mospheres of the Giant Planets (-85%). Therefore it is partic- 
ularly well suited to determine the deuterium abundance from 
the HD/H2 ratio. A number of deuterium abundance determina- 
tions from infrared observations of CH3D have been reported for 
Uranus and Neptune (de Bergh et al.[T986i de Berg h et al il990i 
Orton et al.[T992, Fletcher et al. [2QTo1 Irwin et al. [20T2l) . how- 
ever the knowledge of the required isotopic enrichment factor 
f=(D/H)cH4/(D/H)H2, is uncertain (Lecluse et al. 1996). Direct 
observations of rotational far-infrared HD transitions by ground 
based observatories are difficult because of the opacity of the 
earth's atmosphere in the relevant wavelength regions. Optical 
detections of HD have been reported by Trauger et al. ( 1973 ) and 
Smith et al. (1989a) on Jupiter and Saturn respectively. Similar 
measurements of HD at visible wavelengths, reported for Uranus 
(Trafton & Ramsay "1980) suffer from low signal-to-noise ratio 
and blending with weak CH4 lines, so they can provide only up- 
per limits in the case of both, Uranus and Neptune (Smith et al. 
I1989M . On the other hand, space based observatories like ISO 
(Kessler et al. 11996b and Herschel (Pilbratt et al. [20TQl) have ac- 
cess to the mid- to far-infrared rotational (R-branch) lines of HD 
(Ulivi et al. 1991, see also Table [2]). While the Voyager/IRIS 
instrument could not detect these lines due to insufficient spec- 
tral resolution, observations by the medium-to-high resolution 
spectrometers ISO-S WS (de Graauw et al. [T996l) and ISO-LWS 
(Clegg et al. 11996b resulted in D/H values with significantly 
improved precision for all four Giant Planets (Griffin et al. 
IT996I Encre naz et al. lT996l Feuchtgruber et al. 1999, Lellouch 
et al. '20011 The baseline spectral response calibration of the 
ISO-LWS instrument was obtained from Uranus observations. 
Consequently, this fact prevented a quantitative analysis of the 
HD R(0) 1 12 yum and HD R(l) 57 jim lines for both Uranus and 
Neptune. An analysis of ISO-LWS measurements of Mars and 
Callisto gave only inconsistent results (Davies et al. '2000*), i.e. 
a line detection at low signal-to-noise with a relative spectral re- 
sponse derived from Callisto and a non-detection at high signal- 
to-noise with a relative spectral response derived from Mars. 

The Herschel-PACS instrument is therefore the first astro- 
nomical facility allowing to access the information provided by 
these two lines with high signal-to-noise ratio. 

One major goal of the Herschel Key Program HssO ("Water 
and Related Chemistry in the Solar system", Hartogh et al. 2009 ) 
was indeed to observe HD lines in the four Giant Planet atmo- 
spheres and improve the uncertainties of their respective D/H 
ratios. First results from HD observations on Neptune by the 
far-infrared spectrometer PACS (Poglitsch et al. 2010) on board 
Herschel have been presented in Lellouch et al. (2010). In this 
work we present a combined analysis of the three rotational far- 
infrared lines of HD detected with ISO-SWS (HD R(2)) and 
Herschel-PACS (HD R(0) & HD R(l)) towards a more accurate 
D/H ratio in the atmospheres of Uranus and Neptune. 

2. Observations and Data Reduction 

PACS spectrometer observations of the HD R(0) and R(l) lines 
on Uranus and Neptune have been performed in the high spectral 
sampling density chop-nod mode of the instrument (Poglitsch 
et al. 2010 ). The details of the observations are summarized in 
Table[T] The spectrometer spectral resolution A/SA at the two HD 
line wavelengths of 56.23 yum and 112.07 jim is about 2500 and 



950, respectively. Each of the lines has been measured both in 
range scan mode (covering large wavelength ranges) and in line 
scan mode (short wavelength interval centered on the line). A 
number of further shallow grating scans covering the full PACS 
wavelength range have been carried out on both planets as part 
of the PACS instrument calibration program on their continua, 
however those measurements did not reach sufficient signal-to- 
noise ratios on the two HD lines and are not included here. As 
part of the Herschel Open Time (OT) program "Variability in 
Ice Giant Stratospheres: Implications for Radiative, Chemical 
and Dynamical Processes" led by G. Orton (OTl_gorton01_l 
program), a number of observations of the HD R(0) and HD 
R(l) lines have been taken at diff'erent longitudes. The HD R(l) 
line was part of the observing program on Uranus, but not on 
Neptune. Fortunately this line is seen in the blue spectrometer 
channel simultaneously to the red channel R(0) line observa- 
tions on Neptune, however only in grating order 2, at accord- 
ingly lower spectral resolution of 1400. These data are included 
here as well, allowing us to improve significantly on the statistics 
and to assess observational reproducibility. The data have been 
extracted from the Herschel science archive and processed up to 
Level 1 within HIPE 8.0 (see Poglitsch et al. 2010). Rebinning 
and the combination of the two Nod positions has been done out- 
side HIPE by standard IDL tools. The absolute calibration uncer- 
tainty of PACS spectrometer data is about 30% (Poglitsch et al. 
2010 ). For an accurate determination of the D/H ratio, all spectra 
have been divided by their local continua. Thus, absolute calibra- 
tion errors cancel out and the uncertainties on the line contrast 
are only driven by the signal-to-noise ratio of the observations. 
Although the wavelength scale of the observations has been cor- 
rected by the Herschel-Target velocity Doppler shift, residual 
line center shifts remain. The size of the spatial PACS spectrom- 
eter pixels is about 9.4"x9.4" and the width of the Herschel tele- 
scope PSF (Point Spread Function) ranges from about 6" to 15" 
within the PACS wavelength range. Consequently Uranus and 
Neptune can be considered point sources in the context of these 
observations. The nominal spacecraft pointing uncertainty can 
move point sources at significantly diff'erent positions within the 
PACS spectrometer slit. As a result, the wavelengths of the spec- 
tra may appear slightly shifted, since the nominal calibration ap- 
plies to the slit center or for extended sources only (Poglitsch et 
al. 120 101) . Therefore, to prepare for a best-fit analysis, the spectra 
of all lines have been recentered onto their rest wavelengths. 

Long wards of the R(l) line, around 56.325 yum, a strong 
stratospheric H2O emission line is detected. It is however not 
included in this model, as being inconsequential for the D/H de- 
termination. 

The ISO-SWS observations of the HD R(2) line at 37.7 jum 
have been carried out in 1996-1997 on the two planets at a spec- 
tral resolution of 1700. Together with the HD line, the quadrupo- 
lar rotational lines of H2, S(0) at 28.22 jim and S(l) at 17.03 yum, 
have been measured to obtain independent constraints for the 
thermal profile modeling at similar atmospheric pressure levels. 
Observational details, data reduction and modeling of the HD 
R(2) line on the two planets are described in Feuchtgruber et al. 
(il99 9). To facilitate a common modeling scheme with respect 
to the PACS spectra, these data have also been divided by the 
continuum. 

3. HD Line parameters 

The available information on spectroscopic HD line parameters 
from the literature has been revisited and a few significant up- 
dates with respect to the line parameters used by Feuchtgruber 
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Fig. 1. Monochromatic contribution function for the three HD Hnes (red: R(0); green: R(l); blue: R(2)). SoHd: Hne core; Dashed: 
continuum 



Table 1. Summary of PACS observations 



Target 


Exposure 


HD Line 


Date 


A range 




[sec] 






[fim] 


Uranus 


6941^ 


R(0) 


24-NOV-09 


102-145 


Uranus 


7996^ 


R(l) 


25-NOV-09 


52-62 


Uranus 


2394^ 


R(0), R(l) 


6-Jul-lO 


56, 112 


Uranus 


2x1210^ 


R(0), R(l) 


12- Jan- 11 


56, 112 


Uranus 


2x1210^ 


R(0), R(l) 


13- Jan- 11 


56, 112 


Uranus 


8^x1210^ 


R(0), R(l) 


5-Jun-ll 


56, 112 


Neptune 


6941^ 


R(0) 


30-Oct-09 


102-145 


Neptune 


7996^ 


R(l) 


30-Oct-09 


52-62 


Neptune 


3168^ 


R(0), R(l) 


25-May-lO 


56, 112 


Neptune 


8" x730^ 


R(0) 


5-Jun-ll 


112 


Neptune 


8" x730^ 


R(l) 


5-Jun-ll 


56 



Notes. All observations have been carried out in standard chopped- 
nodded mode. The R(0) line has been measured in the grating or- 
der and the R(l) line in 3^^ one, except the 8 observations of the R(l) 
line on Neptune, where the line was only seen in the 2"^ grating order. 

Range scan mode Line scan mode Observations taken within the 
OTl_gorton01_l observing program have been executed several times 
on the same day at different longitudes to assess longitudinal variability 
of temperature. 



et al. (119991) and Lellouch et al. (120 IQt have been worked out. 
Accurate measurements of the wavelengths of the pure rotational 
transitions of HD are taken from Evenson et al. ( 1988 ) and Ulivi 
et al. (|1991t and remain unchanged. The dipole moment yu of HD 
is taken now from Table II of Lu et al. (119931) as a mean value 



across the 4 measured rotational transitions R(0) to R(3). This 
value of yu=8.21 Debye is then translated into the spectral line 
intensities given in Table [2] A Icr uncertainty of ~3% on the re- 
sulting spectral line intensity values is estimated from the scatter 
in reported dipole moments. Line broadening coefficients y have 
been measured by Lu et al. (199 3 ) at three different temperatures 
compared to 295 K only by Drakopoulos & Tabisz ( |1987al ) and 
Drakopoulos & Tabisz (|1987bl) . The coefficient n describing the 
temperature dependence of the half width at half maximum y 
has been fit to the measured values by: 

r(T) = r(Tref) X (Tref /T)^ with T,,f = 296 K 

Within the temperature range of -50-120 K, containing the 
dominating contributions for HD line modeling in the atmo- 
spheres of Uranus and Neptune, the errors of this fit are <0.5% 
with respect to a linear interpolation of the measurements. For 
the line parameters y, n and the wavenumber shift coefficient S 
of the four transitions, the contributions from the two most sig- 
nificant collisions HD-H2 and HD-He are weighted according 
to the relative abundance of H2 and He in the atmospheres of 
the outer planets (0.85/0.15). Updated rotational constants for 
the HD molecule have been taken from Ulivi et al. ( 1991 ) and 
have been included in the code of Ramanlal &Tennyson ( 20041) 
to calculate the temperature dependence of the partition func- 
tion. The entire set of HD line parameters that has been used in 
the modeling is given in Table |2] The respective values for the 
R(3) transition which is not part of this analysis are provided for 
completeness. 
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Table 2. Adopted HD line parameters for Uranus and Neptune 



Line 


y 


Line intensity^ 






n^ 


6' 






cm~^ / (molec cnT^) 


cm ^/atm 






cm~^ /atm 


R(0) 


89.227950^ 


1.769x10^4 


0.0130 





-0.232 


0.0013 


R(l) 


177.841792^ 


7.517x10^4 


0.0101 


89.228 


0.198 


0.0016 


R(2) 


265.241160^ 


8.870x10^4 


0.0084 


267.070 


0.130 


-0.0045 


R(3) 


350.852950^ 


4.867x10^4 


0.0086 


532.311 


-0.030 


-0.0030 



Notes, n: Exponent of the temperature dependence of the line half width at half maximum y. S: The shift in wavenumber of the line as a function 
of pressure p [atm]. Both parameters are weighted averages according to the relative contributions of HD-H2 and HD-He collisions (~0.85/0.15). 
From Evenson et al. (1988D From Ulivi et al. ( 119911) Derived from Lu et al. ([1993| with T,,/=296K 
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Fig. 2. Thermal (p,T) profiles for both planets. Neptune: dark blue = baseline profile (Bezard et al. 1991); green = Feuchtgruber et 
al. (|1999i) ; light blue = Lellouch et al. (2010); red = best fit profile to the three HD lines; Uranus: blue = baseline profile (Orton et 
al. 120131) : green = Feuchtgruber et al. fl999) ; red = best fit profile to the three HD lines; light blue= best fit to R(0) and R(l) lines 
only 



4. Analysis of the D/H ratio 

The observations were analyzed by means of a multilayer 
radiative-transfer model in which the HD/H2 mixing ratio is as- 
sumed to be uniform with altitude. The model includes atmo- 
spheric opacities of the three measured HD lines (see Table 
|2]) and collision-induced absorption of H2-H2 (Borysow et 
al. 1985), H2-He (Borysow et al. TW), H2-CH4 (Borysow 
& Frommhold 1986) and CH4-CH4 (Borysow & Frommhold 
IT9871) . Orton et al. (2007) updated the H2-H2 coefficients of 
Borysow et al. ([1985 ), but the modifications to far-infrared ab- 
sorption were insignificant. 

The synthetic spectra were calculated monochromatically, 
integrated over all viewing angles of the planets, and then 



convolved with the instrumental profile (FWHM = 120, 215, 
315 km/s) at the respective wavelengths and grating orders. 
Monochromatic contribution functions for the line centers and 
their adjacent continua (Figure [T]) indicate the range of layers 
within the thermal profile dominating the continuum divided 
spectra of the two planets. The line-to-continuum ratios of the 
ensemble of the three HD lines are sensitive to atmospheric pres- 
sure levels from ~10"^ bar to ~1.5 bar. An initial thermal pro- 
file (p,T) of Uranus has been taken from Orton et al. ( 1201 31 ), 
representing a best fit to CH4 and CH3D emission spectra from 
Spitzer IRS data. The initial profile for Neptune is from Bezard 
et al. (119911) , whose tropospheric part relevant for this work orig- 
inates from the Voyager radio occupation experiment (Lindal et 
al. ll990b . These baseline profiles for the two planets do not allow 
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to obtain a good match of all three observed HD lines with the 
model. Together with the HD/H2 mixing ratio, these thermal pro- 
files are adjusted by a 4^^ order polynomial (in log(p)) minimiz- 
ing the least squares to fit all the continuum divided spectra of 
the three HD transitions simultaneously by the model. Figure |2] 
shows the best fit thermal profile for the two planets, together 
with the initial input profiles and earlier work for comparison. 
Combining the Herschel data of the HD lines with Spitzer data 
for a joint modeling is deferred to the future. 

In the case of Uranus, the resulting profile is in remark- 
able agreement with Feuchtgruber et al. (1999 ) and Orton et al. 
(120131 for pressures >200 mbar. The profile of Feuchtgruber et 
al. ([1999 ) is about 5 K colder at 70 mbar, 5 K warmer at 10 mbar 
and about 15 K colder for pressures lower than 1 mbar than our 
best fit profile. With respect to Orton et al. (120131) , the new profile 
is slightly warmer (up to 2 K) between 100 mbar and 5 mbar, for 
lower pressures the new profile becomes continuously warmer 
up to 10 K at 0.1 mbar however. For Neptune, the baseline pro- 
file from Bezard et al. ( 1991 ) and Lindal et al. (1990 ) required 
a slight cooling of 1-2 K starting at 5 mbar towards higher pres- 
sure levels. At lower pressures, the diflferences are then negligi- 
ble. The profile of Lellouch et al. (12010b is warmer at all levels, 
from about 2 K at 1 bar to 4 K at the tropopause up to around 
10 K at 1 mbar. 

In the case of Neptune, the sub-Earth latitude changed only 
from about -27.3° to -28.3° between 1996 and 2010, justifying 
the approach of combining the ISO-SWS data with Herschel- 
PACS data for a joint modeling by one thermal profile. However 
the change of sub-Earth latitude for Uranus went from -46.6° 
in 1996 to -hl3.5° in 2010. VLA mapping observations of 
Uranus at centimeter wavelengths are reported by Hofstadter & 
Butler (2003 ) and Hofstadter et al. ( 201 1 ), which indicate spatial 
changes to the thermal profile at pressures > 1 bar within the pe- 
riod of interest. In order to verify whether a possible change in 
the disk-averaged thermal profile between the epochs of the ISO- 
SWS and PACS observations may aff'ect our analysis on Uranus, 
we derived also a thermal profile by a best fit to the R(0) and R(l) 
lines only. The respective D/H for fitting only these two lines 
is slightly larger (4.7x10"^) but still within the Icr uncertainty, 
however for such a thermal profile the model over-predicts the 
R(2) line already by 30%. The resulting profile (see Fig.|2]) is es- 
sentially the same for pressures larger than lOOmbar compared to 
the profile derived from all three lines. Small departures (<2 K) 
occur only at lower pressures, confirming that our analysis is not 
aff'ected by seasonal changes when including all three lines in 
the calculations. 

Figures [3] and |4] show all the observed spectra and the calcu- 
lated models for Uranus and Neptune respectively. The peak-to- 
peak scatter in line-to-continuum ratios between diff'erent PACS 
observations is about 19% (14 observations) and 17% (10 ob- 
servations) for the R(0) transitions on Uranus and Neptune re- 
spectively. However, this scatter is mainly due to variations in 
spacecraft pointing off'sets, which may cause slight instrumental 
profile variations and therefore variations of the peak contrast. 
A variation of the best-fit D/H ratios by ~ ±1 x 10"^ using the 
same thermal profiles is matching already well to the extremes 
within the R(0) and R(l) observations. However, at the same 
time, variations by this amount appear incompatible with the 
measured spectra of the R(2) lines. The D/H ratio which matches 
the extremes within all observations of the R(0) and R(l) lines 
can be considered as the 3cr statistical error. Translated to Icr 
uncertainties we get D/H=(4.41 ± 0.34) x 10"^ for Uranus and 
D/H=(4.08 ± 0.33) x 10"^ for Neptune. 



However, on top of the statistical error from the number of 
independent observations, there are also systematic uncertain- 
ties on the spectroscopic line parameters of HD. Line intensities 
(3%, Icr), broadening parameter and its temperature coefficient, 
and the wavenumber- shift parameter have independent and tran- 
sition specific uncertainties of a few percent, which may affect 
the calculations of the D/H ratio either way. We therefore add a 
5% uncertainty in quadrature to our statistical error and finally 
quote D/H=(4.4 ± 0.4) x 10"^ and D/H=(4.1 ± 0.4) x 10"^ for 
Uranus and Neptune respectively. Note finally that a little warm- 
ing or cooling of the best fit thermal profiles by +1K leads to 
synthetic spectra which depart already by ~ Icr for all observed 
lines on Uranus and slightly more than Icr for the R(0), R(2) and 
about 2cr for the R(l) measurements on Neptune. Therefore, er- 
ror bars due to thermal profiles uncertainties can be neglected 
compared to those due to the scatter in observed line contrasts. 

These D/H values can be compared to results on the deu- 
terium content in the methane reservoirs of the two planets. 
Using the isotopic enrichment factors from Lecluse et al. (119961) 
of / = 1.68 ± 0.23 and / = 1.61 ± 0.21 for Uranus and Neptune 
respectively with: 

^ (D/H)cH4 
^ (D/H)h, 

and 

CH3D/CH4 = 4 X (D/H)cH4 

we get CH3D/CH4 (Uranus) = (2.96 ^^;^i) x lO'^ and 

CH3D/CH4 (Neptune) = (2.64^^-^^) x 10 -^ in good agree- 
ment with recent results by Irwin et al. C 2012I) (CH3D/CH4 
= (2.9!^;^) X 10-4 for Uranus) and Fletcher et al. (I2QTq1) 
(CH3D/CH4 = (3.0 ^I q) X 10-4 for Neptune). 

5. Discussion 

Our new measurement of the D/H ratio in H2 in Uranus and 
Neptune can be combined with a model of their internal struc- 
ture to constrain the D/H in their protoplanetary ices. Following 
the approach proposed by Lecluse et al. (.1996 ) and also adopted 
by Feuchtgruber et al. (11999b , (D/H)ices can be expressed as: 

(D/H)pianet " XH2(D/H)proto 

(D/H)ices = ^ 

(1 -XH2) 

where (D/H)pianet is the bulk D/H ratio in the planet, taken 
equal to its value in the fluid envelope (D/H)enveiope- This as- 
sumes that the atmosphere and interior of the planet have been 
fully mixed, i.e. that high-temperature equilibration of deuterium 
between hydrogen and ices has occurred during the planet's his- 
tory. This assumption of global mixing is central to our analysis. 
Formation models (Pollack & Bodenheimer 1989 ) suggest that 
the planetary envelopes mixed in early stages, but whether this 
was true also for core material is admittedly unknown. Current 
Giant Planets may not be fully convective, especially Uranus 
(see Podolak et al. 119951 GuiUot 120051) , and it is not known 
whether this state is primordial or not. Another assumption of 
the model is that (D/H)enveiope is equal to the D/H value we 
determine in H2. This assumes that the atmospheric deuterium 
content is largely dominated by H2, with negligible contribution 
from heavy H-bearing species (e.g. H2O). This hypothesis is fur- 
ther discussed below. 

For the protosolar D/H value, we adopt (D/H)pi.oto = 
(2.25±0.35)xlO-^ based on the analysis of ISO-SWS measure- 
ments on Jupiter (Lellouch et al. 2001 ■) . For the volumetric ratio 
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= nH2/(nH2+nH2o) of H2, we initially used the interior mod- 
els of Podolak et al. ( 1995 ). For Neptune, Podolak et al. ( 19950 
considered two model variants, once with a "canonical" density 
in the ice shell, and one with a density reduced by 20 %. These 
models provide values (expressed in Earth masses) for the gas 
(MH2+He). ice (Mice), and rock (Mrock) components of the planet. 
Note that all the Podolak et al. (1995 ) models assumed a so- 
lar rock-to-ice ratio (~2.5), i.e. an ice mass ratio in the heavy 
element component (F = Mice/(Mice + Mrock)) equal to 0.715. 
Assuming solar composition , i.e. H 2/(H2+He) = 0.747 by mass 
(consistent with Helled et al. 1201 lb , xh2 can be expressed as: 

1 



(mH20/'WH2)X/H2) 

where 

0.747M^,^^, 

^""^ 0.141 Mn.^He + Mice 

is the mass ratio of H2 and mH20 and mH2 are the mo- 
lar masses of H2O and H2 respectively (i.e. 18g and 2g). The 
XH2 values and the inferred (D/H)ices are given in Table |^ (note 
that the xh2 values are slightly different from those given in 
Feuchtgruber et al. (11999 ), due to a minor mistake in that pa- 
per). 

We also used the more recent Uranus and Neptune models of 
Helled et al. (120111) . These models propose empirical pressure- 
density models tuned to match the planetary radii, masses, grav- 
itational coefficients J2 and J4 and solid rotation periods deter- 
mined by Voyager. The models are then interpreted in terms of 
bulk composition, i.e. the mass fractions of hydrogen (X), he- 
lium (Y), and heavy elements (Z). Two limiting cases are con- 
sidered for heavy elements, being either pure rock (represented 
for definiteness by Si02) or pure ice (represented by H2O). 
Furthermore, two variants (I and II) are considered regarding 
the radial distribution of the heavy elements within the plane- 
tary interior. Obviously the pure Si02 cases are excluded from 
the point of view of the D/H ratio, since they would lead to a 
protosolar D/H. We thus considered here the H2O cases. By def- 
inition, these models have F = 1 . In this case, fH2 is simply given 
asfH2 =X/(X-fZ). 

Even more recently, interior models of Uranus and Neptune 
were updated by Nettelmann et al. (2013 ). These authors provide 
results based on full sets of three-layer interior models combined 
with different solid-body rotation periods, gravitational data and 
physical equations of state (2 models for Uranus and 3 models 
for Neptune). Values for xh2, fH2. and resulting (D/H)ices from all 
their models are given in Table [3] 

For all the interior models considered. Table [3] indicates 
that the D/H values for the protoplanetary ices of Uranus and 
Neptune are consistently 4-6 times smaller than those found in 
the water-ice reservoirs of Oort-cloud comets (~(2-3)xl0"^, see 
Bockelee-Morvan et al. 2012), and still a factor of 2-3 lower 
compared to the D/H in Earth's oceanic water (1.5x10""^), car- 
bonaceous chondrites (1.4x10"^) and the Jupiter- family comet 
103P/Hartley 2 (1.6x10-4, Hartogh et al.'20lT). 

It is somewhat unexpected to find that proto-uranian and 
proto-neptunian ices are much less D-rich than cometary ices, 
warranting a brief discussion of the above hypothesis and deriva- 
tion. A first issue is that (D/H)enveiope may not be equal to 
(D/H)h2. This situation occurs if the planet envelope is heav- 
ily enriched in D/H-rich volatiles. However, even for an extreme 
water enrichment within the atmosphere (e.g. 0/H = 440 times 
solar, as advocated for Neptune by Lodders & Fegley (.1994 J ), 



Table 3. Inferred D/H in protoplanetary ices 



Model 


fH2 


XH2 


(D/H)ices 


Uranus^ 


0.108 


0.521 


(6.75+1. l/-1.2)xl0-5 


Neptune V 


0.133 


0.581 


(6.7+1.4)xl0-^ 


Neptune 2^ 


0.055 


0.343 


(5.1+0.75)xlO-5 


Uranus i (H2U) 


U.Uo / 


\JAo3 


(d.Zj + 1.Uj)X1U 


TTriiniic TT 

UldllUIS ii \n2\Jj 


U.UD J 


U.JO / 




Neptune I (HsO)^ 


0.082 


0.444 


(5.58+l)xlO-^ 


Neptune II (HsO)^ 


0.074 


0.417 


(5.42+0.95)xl0-5 


Uranus V 


0.116 


0.541 


(6.94+1.3)xlO-5 


Uranus 2^ 


0.097 


0.492 


(6.48+1. 1)x10-5 


Neptune V 


0.150 


0.614 


(7.04+ 1.6)x 10-5 


Neptune 2a^ 


0.178 


0.661 


(7.71+1.9)xl0-5 


Neptune 2b^ 


0.124 


0.560 


(6.46+ 1.4)x 10-5 



Models from: (a) Pod olak e t al. ([T9951) . (b) Helled et al. ([20TT1) 
(c) Nettelmann et al. ([20T3l) 
Volumetric ratio: = nH2/(nH2+nH2o) 
Mass ratio: = Mh2/(Mh2 + Mice) 



the (D/H)enveiope is increased only by a factor 1.24 ((Lecluse et 
al. 1996), leading to only a -30 % increase of (D/H)ices, by far 
insufficient for reconciliation with cometary values. The second 
possibility is that the assumption of complete mixing is not valid. 
In a state of incomplete mixing, the derived values of (D/H)ices 
as per Table [3] would represent lower limits of the D/H in the icy 
cores. Note however that the indistinguishable values of (D/H) in 
H2 (within error bars) in Uranus and Neptune (which still differ 
radically in terms of their internal energy sources and therefore 
presumably of their convective state) suggest that the current in- 
complete convection in Uranus' interior may not be important 
in this respect. Therefore, we feel that this scenario is also not 
obviously promising. 

The (D/H)ices values we infer for the protoplanetary ices cor- 
respond to a modest enrichment factor f~2-3 over the proto- 
solar value. Compared to the most pristine (i.e. D-rich) solar 
system material (i.e. the D/H-rich component of the LL3 me- 
teorites, having f=35) or even to typical cometary material (f~7- 
20), this implies that the protoplanetary ices have been consid- 
erably reprocessed in the solar nebula. Yet evolutionary mod- 
els (Kavelaars et al. 120111 ) accounting for radial turbulent mix- 
ing within the nebula predict much larger enrichments (f=14- 
20) at the estimated 12-15 AU formation distance of Uranus 
and Neptune. Along with the recently revealed diversity of the 
D/H ratio in comets (Hartogh et al. 2011, Bockelee-Morvan et 
al. 120 12t and the apparent absence of correlation of the values 
with the estimated formation distance of these comets (i.e., the 
Kuiper-Belt vs Oort cloud families) our result of a low D/H in 
the Uranus and Neptune original ices would illustrate the limita- 
tion of these evolutionary models. Note also that given the error 
bars, we are unable to find any significant difference between the 
D/H in the proto-uranian and proto-neptunian ices (see Table [3]), 
so searching for a correlation between the D/H and the formation 
distance (as predicted by the evolutionary models), is precluded. 

The idea that the protoplanetary ices should necessarily have 
a deuterium content equal to that measured in cometary water 
may however be challenged by the work of Alexander et al. 
(2012). They find a linear correlation between the D/H and C/H 
measured in a set of carbonaceous chondrites (CC) from the CM 
and CR groups that experienced different degrees of aqueous al- 
teration. This suggests that the hydrogen isotopic composition 
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in these bodies results from the mixing between hydrated sili- 
cates and organic matter, and extrapolating the relationship to 
C/H = should therefore give the isotopic composition of wa- 
ter. This approach provides (D/H)h20 -9x10"^ for CMs (but 
-17x10"^ for CRs). In a second approach, the authors correct 
the D/H measured in other types of chondrites from the con- 
tribution of organic material and infer that in addition to CMs, 
the water compositions of the CIs, CO, CV and Taglish Lake 
meteorites are less deuterium-rich than comets, with D/H gener- 
ally below ~ 10x10"^. As chondrites are fragments of main-belt 
asteroids (and in particular, CCs are generally associated with 
C-type asteroids, which may have formed in the same region as 
comets (Walsh et al. 2011) ), this is thus indication that there 
may have existed material with "low" D/H in water ice origi- 
nating from the formation region of comets. It might therefore 
not be irrelevant to relate the low values we infer for the proto- 
uranian/proto-neptunian ices ((4-9)xl0"^, see Table (sj to such 
material. However, a complication in this scenario is to under- 
stand why so far no comets have been observed to exhibit such 
a low D/H, and what is the origin of the cometary D/H enhance- 
ment compared to this value. 

Coming back to the scenario where cometary ices are repre- 
sentative of protoplanetary ices, all the above calculations rely 
on interior models of Uranus and Neptune, which are not well 
constrained (see e.g. a discussion on the effect of the uncertainty 
of the rotation period by Podolak & Helled 2012 ). It is therefore 
worthwhile to "invert" the problem, i.e. assume some value of 
(D/H)ices and constrain the interior structures. As the pure Si02 
models of Helled et al. (12011 1) cannot explain a (D/H)h2 larger 
than the protosolar value, while the pure H2O models lead to too 
small values for (D/H)ices, it is clear that intermediate models 
(i.e. a mix of ice and silicates) are needed. We therefore searched 
for the ice mass ratio in the heavy element component (F = 
Mice/(Mice + Mrock), targeting (D/H)ices = 1.5x10-4 or SxlO'l 
For a given input value of F, the X, Y, Z values were interpolated 
from Table 3 of Helled et al. (2 011) , considering the average of 
cases I and II. Zice and Zj-ock are then given by FxZ and (l-F)xZ, 
respectively, and (D/H)ices was obtained as before, using now f//2 
= X/(X-hZice). Finally, the knowledge of X, Zice and Zrock permits 
us to derive the mass of H, the total mass of O, and the mass of 
O contained in the ice, which can be translated into 0/H ratios. 
Results are given in Table |4] where the 0/H ratios are expressed 
in mass, volume, and in the enhanced factor over the solar value 
(assuming volume solar 0/H = 4.9x10""^; Asplund et al. ( 2009 )). 
Of course the precise values given in Table |4] are somewhat de- 
pendent on the simplified description, attached to the Helled at 
al. 1201 II models, that all the ice is in the form of H2O and all 
the rock in the form of Si02. And again, results are subject to 
the validity of the complete mixing hypothesis: for incomplete 
mixing, the F values reported in Table |4] would represent lower 
limits to the actual ice mass fraction. 

The models in Table |4] have F = 0.14-0.32. In other words, 
we infer that 68-86 % of the heavy component consists of rock 
and 14-32 % is made of ice. Therefore, unlike in the Podolak 
et al. ([1995 ) models and in the ice (H2O) version of the Helled 
et al. (I2OII ) models, we find that Uranus and Neptune interi- 
ors might be more rocky than icy. This behaviour is similar to 
the case of Pluto, which based on the body density (~2 g cm"^), 
has an estimated rock mass fraction of about 0.7 (Simonelli & 
Reynolds [1989]). All the above models have Z = 78-81 % per 
mass. With Uranus' and Neptune's masses equal to 14.5 and 17.1 
Earth masses, this gives ~1 1.6 of heavy elements for Uranus 
and -13.5 M® for Neptune. This is to be compared with the 
findings by Owen & Encrenaz (120031 12006t who similarly used 



Table 4. "Inferred" interior models 
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3.08 
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2.91 


2.58 


Volume 


0.193 


0.169 


0.182 


0.161 


X solar"" 


393 


345 


372 


329 


0/H ratio (0 in ice) 










Mass 


1.35 


0.619 


1.160 


0.532 


Volume 


0.084 


0.039 


0.073 


0.033 


X solar"" 


172 


79 


148 


68 



Notes. solar O/H volume ratio = 4.9xl0-\ Asplund et al. ( [20091 



the enrichment in heavy elements in the Giant Planets (mostly 
measured in carbon) to estimate the mass of the "SCIP" (solar- 
composition icy planetesimals) within each planet. For Neptune, 
Owen & Encrenaz (2006) find a SCIP mass of (13+3) M®, fully 
consistent with our value. The agreement is worse at Uranus, 
where Owen & Encrenaz (120061 find (8.5+|^) M®. 

The total 0/H ratio is 329-393 times solar, but when only 
O from H2O is considered, the 0/H ratio is 68-172 times solar. 
Independent constraints on the atmospheric 0/H ratio have been 
inferred from the measured CO mixing ratio, a disequilibrium 
species whose tropospheric abundance is sensitive both to the 
vigor of vertical mixing from the deep atmosphere and to the 
0/H ratio (CO is produced from H2O from the net reaction CH4 
+ H2O ^ CO 3H2.) To explain the ~1 ppm CO abundance 
initially measured in Neptune's troposphere (Marten et al. . 19931 
Guilloteau et al. fT993l) Lodders & Fegley ([19941) invoked a 440 
times solar 0/H ratio in Neptune's deep atmosphere, using a so- 
lar 0/H = 7.4x10-4. Rescaling this to 0/H = 4.9x10"^ gives 0/H 
= 660 times solar, which is 4.5-9 times higher than we infer. We 
conclude that the Lodders & Fegley (1994 ) models are incon- 
sistent with our D/H measurement, although reconciliation may 
be possible if the CO abundance is actually overestimated by a 
considerable fac tor. A s a matter of fact, subsequent observations 
(Lell ouch e t al. [20051 Hesman et al. [20071 Luszcz-Cook & de 
Pater 120131) all indicate that Neptune CO has two distinct com- 
ponents and that its tropospheric abundance is lower than pre- 
viously thought. Luszcz-Cook & de Pater (120131) in particular, 
determined a much smaller (-0.1 ppm) CO tropospheric abun- 
dance. However, by reassessing the Lodders and Fegley ([19941) 
model, especially in terms of (i) the limiting reaction steps and 
(ii) the characteristic mixing time, they still find that that a global 
0/H enrichment of at least 400, and likely more than 650, times 
the protosolar value is required to explain their measured CO 
abundance. Therefore, the discrepancy with our estimate of the 
atmospheric 0/H remains. 

Finally, we note that the 0/H volume ratios for the H2O 
component are in the range 0.033-0.084, i.e. H2O / H2 = 0.07 
- 0.17 in the atmosphere. This induces only a minor correc- 
tion to (D/H)h2 (Lecluse et al. 1996), i.e. (D/H)enveiope = (1-04- 
1.09)x(D/H)h2 and does not impact any of the above conclu- 
sions. 



7 



Feuchtgruber et al.: D/H in Uranus and Neptune 



Uranus 



Neptune 
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37.80 



Fig. 3. Observed and synthetic Uranus spectra. Black: contin- 
uum divided spectra (solid: range scan; dashed: HssO line scan; 
dotted: OT line scans) and best fit model (red) with D/H=4.41 x 
10"^ . Synthetic spectra for diff'erent D/H ratios but for the same 
thermal profile to illustrate the sensitivity to this model param- 
eter are shown as green solid lines for D/H=3.4 x 10"^ and 
D/H=5.4 X 10-5 (3cr) and blue solid lines for D/H=4.07 x lO'^ 
and D/H=4.75 x 10"^ (Icr). The spectral line around 56.33 jim is 
caused by stratospheric H2O emission and is not included in the 
model. The R(2) line has been measured by ISO-SWS. 
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Fig. 4. Observed and synthetic Neptune spectra. Black: contin- 
uum divided spectra (solid: range scan; dashed: GT line scan; 
dotted: OT line scans) and best fit model (red) with D/H=4.08 x 
10"^ . Synthetic spectra for diff'erent D/H ratios but for the same 
thermal profile to illustrate the sensitivity to this model param- 
eter are shown as green solid lines for D/H=3.09 x 10"^ and 
D/H=5.07 X 10-5 Oct) and blue soHd lines for D/H=3.75 x lO'^ 
and D/H=4.41 x 10"^ (Icr). The R(l) line has been observed in 
two diff'erent grating orders (n) and accordingly at diff'erent in- 
strumental resolution. The strong spectral line around 56.33yum 
is caused by stratospheric H2O emission and is not included in 
the model. The R(2) line has been measured by ISO-SWS. 



6. Conclusions 

Herschel-PACS spectrometer observations combined with pre- 
vious ISO-SWS observations of the three lowest energy ro- 
tational lines of HD on Uranus and Neptune have been an- 
alyzed. Applying latest spectroscopic line parameters for the 
HD R(0), R(l) and R(2) lines in multilayer atmospheric radia- 
tive transfer calculations, a least-squares fit of synthetic spec- 
tra to the continuum divided measurements results in D/H val- 
ues of (4.4±0.4)xl0-5 and (4.1±0.4)xl0-5 (Icr) for Uranus and 



Neptune respectively. The simultaneous modeling of the three 
measured HD lines requires only small modifications to thermal 
profiles (p,T) known from earlier work on both planets. The new 
D/H values, although somewhat smaller than in previous analy- 
sis, confirm the enrichment of deuterium in the atmospheres of 
Uranus and Neptune compared to the protosolar nebula. Based 
on published interior models in which the largest fraction of the 
heavy elements is in the form of ice (70-100 %), and assum- 
ing that complete mixing of the atmosphere and interior occured 
during the planets history, the required D/H in the protoplanetary 
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ices responsible for this enrichment is significantly lower than 
known from any water/ice source in the solar system. A possible 
solution to this unexpected result is that the interiors of Uranus 
and Neptune are actually dominated by rock. For icy material in 
the form of H2O and rocky material in the form of Si02, an ice 
mass fraction of only 14-32 % is found. Still, a complete inter- 
pretation of the now accurate values of Uranus and Neptune D/H 
would greatly be aided by a more definite picture of their inter- 
nal structures. The latter would benefit from improved gravity, 
shape and rotation data, that could be obtained from planetary 
orbiters. 
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